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Abstract 

We show that the solution pubhshed in Ref.l is geodesically com- 
plete and singularity-free. We also prove that the solution satisfies 
the stronger energy and causality conditions, such as global hyper- 
bolicity, causal symmetry and causal stability. A detailed discussion 
about which assumptions in the singularity theorems are not fulfilled 
is performed, and we show explicitly tat the solution is in accordance 
with those theorems. A brief discussion of the results is given. 

PACS: 04.20.Jb, 04.20.Cv,98.80.Dr 
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1 Introduction 



Recently, one of us presented a solution to Einstein's field equations for 
a perfect-fluid energy-momentum tensor without any curvature singular- 
ity (Ref.l). The solution has positive pressure and energy density every- 
where, and it satisfies a realistic equation of state for hot radiation-dominated 
epochs. As we shall see, it also verifies the stronger causality requirements. 
Therefore, the question remains of whether or not the solution is geodesically 
complete and, if it is, how that fits in with the general conclusions of the very 
powerful singularity theorems [2,3]. 

The main result in this paper is that the mentioned solution is geodesically 
complete and, in fact, singularity-free (see Ref.2 for definitions). We shall 
prove it in Sec. 2 by studying the equations of the non-spacelike geodesies and 
showing that they can always be extended to arbitrary values of the affine 
parameter. Once this important property has been established, in Sec. 3 we 
shall study which of the assumptions of the different singularity theorems are 
not verified in the solution. Since the energy and causality conditions will be 
proven to be satisfied, it will turn out that the so-called initial or boundary 
condition is the one that fails to hold. This condition usually refers to the 
existence of a causally trapped set [2,3] or to some fixed bound for the initial 
expansion of the geodesies. Thus, we shall see that the initial condition is 
essential for the development of singularities. Along with this analysis, we 
shall prove a number of nice and desirable properties of the solution, among 
them, global hyperbolicity and stable causality. Finally, we devote Sec.4 to 
give a brief discussion of the results herein presented. 

Before getting into the main subject, let us give a summary of the solution 
of [1]. The line element is given by 



ds^ = C\at)C\3ar){-dt^ + dr^) + {Qa^)-^C\at)S^{?,ar)C-^'\'iar)d(t)''+ 



where a 7^ is an arbitrary constant (which can be taken as positive) and 
we use the notation 




(1) 



S{u) 



sinh{u), C{u) 




The coordinates used are of cylindrical type, and their range is 



—00 < t, z < 00, 0<r<oo, < (j) < 2% 
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such that and 0+27r are identified. The metric (1) is a solution of Einstein's 
equations for a perfect fiuid. The coordinates are comoving (or adapted to 
the fluid congruence) in the sense that the velocity vector has the form 

u = C-\at)C-\3ar)d/dt. (2) 

The fiuid congruence is orthogonal to the spacelike hypcrsurfaccs t = const. 
and therefore its rotation tensor vanishes. The expansion of the fluid has the 
expression 

e = 3aS{at)C-%at)C-\3ar), (3) 

which is regular everywhere and positive for t > (expanding phase) or 
negative for t < (contracting phase). We also see that both for t — > ±oo 
or r — > oo the expansion approaches zero. 

However, the congruence of the fluid is not geodesic, the acceleration 
vector fleld being 

a = 3aS{3ar)C-^{3ar)C-\at)d/dr. (4) 

From a physical point of view, this fact is of extreme importance for the 
absence of singularities in the solution. Because, as is well known, a non- 
vanishing acceleration means that there exists a spacelike gradient of pressure 
or equivalcntly, a force which opposes to the gravitational attraction. Thus, 
although for t < there is a contraction, the gradient of pressure docs not 
allow for the formation of a singularity and causes a bounce from which the 
matter expands. The fluid congruence is also shearing, the shear tensor can 
be found in [1]. 

The density and pressure of the fluid take the following form 

XP = 15a^C-\at)C-\3ar), (5) 
P = P/3 (6) 

where x is the gravitational constant. From (5) we see that the maximum 
value of p at each time t is achieved at r = 0. Moreover, the absolute 
maximum of the density (and pressure) occurs when r — and t — 0, and 
its value is pm — 15a^/x. Therefore, we can interpret the constant a as 
the maximum energy density. Note that, since a is arbitrary, the maximum 
density can be chosen as large or small as we like. On the other hand, from 
(6) we learn that the equation of state is realistic for a space-time filled 
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with incoherent radiation Equations (5) and (6) also tell us that p and p are 
positive everywhere, which implies that the strong energy condition (or the 
timelike convergence condition, see [2]) is satisfied. Prom the same equations 
we can check that p and p are finite and well-behaved everywhere, which 
means that there is no matter singularity (singularity in the Ricci tensor). 
It can further be shown, by compu ting the Wcyl tensor (see Rcf.l), that all 
the curvature invariants are regular on the whole space-time so that there is 
no curvature singularity in the metric at all. 

The line element (1) admits an abelian group of symmetries G2, and 
the Killing vectors d/dcf) and d/dz are globally defined. Both of them are 
spacelike, orthogonal to each other and orthogonally transitive, and therefore 
the metric is cylindrically symmetric. The regularity condition [4] on the axis 
r = is satisfied so that the 27r-periodicity of the coordinate is well defined 
and the so-called elementary flatness on the vicinity of the axis is assured. 
In this sense, the coordinate singularity which appears in (1) at r = is just 
that of cylindrical coordinates; other coordinates (Cartesian-like) do exist 
such that they provide a global chart on the manifold. 

The coordinate t appearing in the metric is a time function (or cosmic 
time) in the sense that it (decreases) increases along every (past-) future- 
directed non-spacclikc curve. This can be seen by checking that the gradient 
of t is timelike everywhere. As is well known, this is the necessary and 
sufficient condition for the stable causality condition to hold [2]. Thus, the 
space-time is causally stable which is the stronger causality requirement and 
it implies the weaker chronology and causality conditions. 

Finally, let us remark that the solution is inextendible, which avoids the 
possibility of the existence of singularities bidding somewhere. This is ap- 
parent from the form of the metric (1); however, the best way to prove it is 
to show that all geodesies are complete and then the metric can be extended 
nowhere. In the next section we shall prove geodesic completeness for the 
solution. 

2 Geodesic completeness 

Our aim now is to show that every non-spacelike geodesic can be extended 
to arbitrary values of its affine parameter. To that end, we must analyze 
the equations of the geodesies in the metric (1), which after standard and 
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straightforward calculations can be written as follows 



i+2aT{at){P + r^) + 6aT{3ar)ir + 
(2/9a)T{at)S\3ar)C-^/^(3ar)^^ - 

aC-\at)S{at)C-^/^{3ar)z^ = 0, (7) 

r + 3aT(3ar)(i^ + P) + AaT{af)if - 
{l/9a)S{3ar)C-^l^{3ar)[3 - T^{3ar)]^'^ + 

aC-^{at)S{3ar)C-^^l^{3ar)z^ = 0, (8) 

{9a^)-^C\at)S^{3ar)C-^/\3ar)4) = K, (9) 
C-2(ai)C-2/'(3ar)i = L, (10) 

C\at)C\3ar){-t' + f^) + 
(9a2)-iC4(at)^2(3ar)C-2/3(3ar)02 + 

C-\at)C-^'\3ar)z^ = -5. (11) 

Here a dot means derivative with respect to the affine parameter and we use 
the notation 

T{u) = tanh{u). 

K and L are constants of motion along the geodesies due to the existence of 
two Killing vectors in the spacetime. Another constant of motion is denoted 
by 5, which must be taken equal to one or zero for timelike or null geodesies 
respectively. The fact that all functions involved in eqs.(7)-(ll) are non 
singular will grant us the existence and uniqueness of solutions. 

The method of the proof is to give finite bounds for the first derivatives 
of the coordinates, which implies (see Ref.5) that the field is non-singular 
and the geodesies are complete. We will also use the second derivatives of 
the coordinates to show that these cannot become singular. We shall only 
deal with geodesies propagating towards the future. Propagation towards 
the past can be treated in a similar way. 

In equation (8), when (9) and (10) are taken into account, f is negative 
for positive t and increasing values of r. This follows from the fact that the 
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positive term in dominates over the negative term in 0^ for large enough 
values of r. Negative values of t need not be considered, as t will end up 
becoming positive. It is obvious that r cannot diverge to infinity in a finite 
proper time with negative r. 

Most of the geodesies approach r = without any problem. As r de- 
creases, the term in tr dominates over the terms in and t^, while the term 
in vanishes when approaching r = 0. Hence, f is positive for decreasing r 
in the vicinity of r = for positive t. This means that the geodesies cannot 
come near the axis quickly enough to become singular. 

Also t is negative for large values of t (notice that the term in tr would 
not be relevant if r decreases), so that the geodesic does not grow faster than 
its tangent. The quantity z does not become singular if r and t do not; 
could have problems if r were zero, but we shall see that this is not possible. 
In order to proceed, we divide our study in steps, starting with the simpler 
geodesies and then going to the general case. In that way we shall be able 
to use some results several times as we progress on. 

1) Geodesies in the fluid eongruenee. This is the simplest case. The 
congruence of the fluid is defined hy f — (f) — z — Prom Eq. (4) it follows 
that the only possible geodesies he in r = 0. Then the system reduces to the 
single equation 

i = C-^{at) < 1. 

Obviously, these geodesies are complete. 

2) Geodesies along the axis. They are defined byf = = Or = O. The 
remaining equations are now 

z = LC^{at), 

i = C-\at)[L^C\at) + 5]'^' < (L^ + Sf/^. 

Therefore, the derivatives are bounded everywhere and these geodesies are 
also complete. The coordinate t reaches infinity only when so does the affine 
parameter. This is a general property of all geodesies and we shall not refer 
to it in what follows. 

3) Radial null geodesies. Now (f) — z — 0, 5 — 0. After one integration 
the system reduces to 

i=\r\, r = hC-\at)C-^{3ar) 
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where h is a constant of integration. We see that < |r| = \t\ < \h\, 
from where it follows that these geodesies are complete too. When one of 
the curves arrives at r = along a direction (p then it continues along the 
direction + tt after crossing the axis without problems. 

4) Radial timelike geodesies. Again = i = 0, but now S = 1. By 
introducing a new function v which, so to speak, parametrizes t and r, the 
system reduces in this case to 

i^C{v)A{r,t), S{v)A{r,t), v ^ B{r,t,v)A{r,t) (12) 

where 

A{r,t) = C~\at)C-\3ar), (13) 

B{r, t, v) = -a[3T{3ar)C{v) + 2T{at)S{v)]. (14) 

The previous reasoning about second derivatives allows us to conclude that 
these geodesies are complete. 

5) Null geodesies with no angular velocity. Now, = and 5 — 0. Using 
the function v again, we can write the system of equations 

i^C{v)E{r,t), r ^ S{v)E{r,t), v ^ E{r,t)F{r,t,v) 

where 

E{r,t) = \L\C-\at)C-'^/^{3ar), (15) 
F(r, t, v) = -a[4:T(3ar)C(v) + 3T{at)S(v)] . (16) 

A reasoning similar to that given in the previous case leads us again to 
completeness for these null geodesies. 

6) Null geodesies on the hypersurfaces z— const.. These are defined by 
z — 0, 5 — 0. The system of equations can be written as 

i = C{v)M{r, t), r = S{v)M{r, t), v = M{r, t)D{r, v) 

with 

M(r,t) = 3\aK\C-\at)C-^/^{3ar)S-\3ar), (17) 
D{r,v) = aC{v)[3T-\3ar) - 4:T{3ar)]. (18) 
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In this case we can obtain one of the orbit equations by dividing r by v and 
then integrating. Doing so we get 

C{v) = a-^S{3ar)C-'^/\3ar) 

where a is a positive constant. We see that, since C{v) > 1, the co- 
ordinate r can only take values between r+ and r_, with r± defined by 
S'(3ar±)C~^/^(3ar±) = a. Thus |v| is bounded and the geodesies are com- 
plete. 

Among these geodesies there is one which is circular and goes round the 
axis of symmetry at r = = arccosh{2)/3a. This can be identified as the 
special case in which r+ — — ro, v — 0. 

7) General non-spacelike geodesies. Using the f-parametrization, we can 
rewrite the general system of equations (7)- (11) in the following form 

i = {3\aK\)-^M{r, t)G{r, t)C{v), r = {3\aK\)-^M{r, t)G{r, t)S{v), 

V = {3\aK\)-^M{r, t)H{r, t)G-\r, t) (19) 
where we have defined 

G{r,t) = [9a^K^ + L''C%at)S\3ar) + C\at)S\3ar)C-^/''{3ar){l - 6)f\ 

(20) 

H{r,t,v) = -a{S{v)T{at)S\3ar)[3L^C%at) +2C\at)C-^/\3ar){l - 6)] + 
C{v)T{3ar)[4L^C%at)S\3ar) + 3C\at)S^{3ar)C-^/\3ar){l - 6)] + 

9a^K^C{v)[4:T{3ar) - 3T^\3ar)]}. 

(21) 

The only problem wc could have appears when r approaches zero, but in 
this case the dominant terms are those of paragraph 6), so that the geodesies 
with K ^ cannot reach the axis. Then they are complete. We can say that 
the term in 0^ in (8) is centrifugal in the vicinity of the axis. 
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3 Properties of the solution and its relation 
with singularity theorems 

Prom the results of the previous section it is obvious that every maximally 
extended null geodesic meets any of the hypersurfaces t — const.. This 
is a sufficient condition such that every non-spacelike curve intersects the 
mentioned hypersurfaces exactly once [6]; in other words, the hypersurfaces 
t = const, are global Cauchy surfaces and the solution is globally hyperbolic. 
For this reason the solution is also causally simple [2], that is, for every 
compact set K, J'^{K) is closed and its boundary coincides with E'^{K). 

Due to the fact that every spacelike hypersurface t = const, is a global 
Cauchy surface, and since t is a time function in the solution as explained 
in Sec.l, it follows that every non-spacelike curve (geodesic or not) can be 
extended to arbitrary values of its generalized affine parameter (because it 
has to meet all the Cauchy surfaces). This means that the solution is non- 
spacelike 6-complete (see Ref.2), that is to say, singularity-free. 

Furthermore, from expressions (5) and (6) we see that the energy-momentum 
tensor does not vanish anywhere. Therefore, given any non-spacelike vector 
V we have RabV^'v'' > so that, in addition to the strong energy condition, the 
generic condition is fulfilled as well [3]. Thus, it remains to see which other 
conditions in the singularity theorems arc not satisfied by the spacetime. 

The simplest such theorem (Ref.2, p. 274) requires the existence of a 
Cauchy surface such that the trace of its second fundamental form is bounded 
away from zero. In our case, the mentioned trace for any Cauchy surface 
t — const, is given by the expansion of the fiuid congruence shown in (3). 
From this expression is obvious that is not bounded away from zero, be- 
cause, for any t, 9 goes to zero when r — oo. The significance of this result 
can be better understood as follows: given any Cauchy surface t = ti < 0, we 
have 6t^ < and then there are points conjugate to the surface along every 
future-directed timelike geodesic orthogonal to the surface within a distance 
— 3/^tj. However, as 6t^ is not bounded away from zero, those distances do 
not have an upper bound. This allows for the existence of a maximal geodesic 
orthogonal to the surface up to any point q to its future, as is necessary being 
t — ti a, Cauchy surface. 

Another property of the solution is that the Cauchy surface E : i = is 
a maximal spacelike hypersurface, i.e., one in which the trace of the second 
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fundamental form vanishes. It follows that the spacetime is causally (or time) 
symmetric [7] so that J+(S), /^(S) and -D^(S) are isometric to J^(S), /^(S) 
and -D~(E) respectively. In Ref.7, Tipler has shown that all geodesies are 
both future and past incomplete in a spacetime which contains a non-compact 
maximal Cauchy surface if, in addition to the strong energy condition, the 
following is satisfied: there exist fixed positive constants b, c such that 



for every timelike geodesic (with tangent vector v) intersecting E orthogo- 
nally at r = 0, where r is the afhne parameter. However, this is not verified 
in solution (1), because given any pair of constants b and c, the above integral 
is always positive but not bounded below above zero. In fact, we can always 
choose geodesies (for r big enough initially) such that the integral takes a 
value as small as we like and less than any previously fixed c. 

Let us consider then the classical singularity theorems. The first of them, 
chronologically speaking, is that of Penrose [8,2], which assumes the existence 
of both a non-compact Cauchy surface and a closed trapped surface, i.e. a 
compact (without boundary) spacelike two-surface in which the traces of the 
two null second fundamental forms have the same sign. Our space-time does 
not have a closed trapped surface. In order to prove it, suppose there was 
one. Since the surface is compact, it must have a point p where r reaches its 
maximum, so that the normal at p is a superposition of d/dr and d/dt at p. 
But if we compute the traces of both null second fundamental forms at p, we 
get 



since n^^zip) and n^,p{p) are positive for outgoing normals and negative for 
ingoing normals at p. We see that the traces have opposite signs so that 




2-^/^aC-^(at)S-\3ar)[-T{at)T{3ar) T 3 ± 2T\3ar)] 



from where it is clear that 



X-d< -2'^^^aC-^{at)S-\3ar)[l + T{at)T{3ar)] < 0, 
xid > 2-^/^aC-\at)S-\3ar)[l - T{at)T{3ar)] > 0, 



10 



there are no closed trapped surfaces. This reasoning can be visuahzed in the 
following way. Take any closed compact two-surface in the manifold. This 
surface must be orthogonal to the in- and out-going radial null geodesies 
somewhere. But the outgoing and ingoing radial null geodesies are expanding 
and contracting, respectively, everywhere. 

The most famous and powerful singularity theorem was proven by Hawk- 
ing and Penrose [9,2]. Our space-time satisfies all conditions in the Hawking- 
Penrose theorem except for condition 4 as appears in Ref.2, which allows for 
three different possibilities. One of them is the existence of a closed trapped 
surface, which we have just shown is not satisfied in the solution. The second 
possibility is the existence of a point q such that on every past (future) null 
geodesic from q the expansion becomes negative. There is no such point in 
solution (1), since through any point in the manifold there are radial null 
geodesies which diverge if they are outgoing and future-directed or ingoing 
and past-directed. Another way to see that the point q does not exist is 
to remember that through any q there are null geodesies with z = const. 
which are bounded above and below in r. Thus, these geodesies can never 
converge with the radial ones. The third possibility in Hawking-Penrose's 
theorem is the existence of a compact achronal set without edge. It is pretty 
obvious that there is not a set with those properties in the space-time under 
consideration. To prove it rigourosly, take any achronal set in the mani- 
fold and a point q in the set. By using the radial geodesies we can always 
choose points q- e I^iq) and q^ e ^'^(q) such that r(g_) = r{q+) > r{q), 
where for any point s we denote by r(s) the value of the coordinate r at 
s. Since q^ G and r(g+) = r{q_), we can join g_ and g+ with a 

future-directed worldline of the fluid congruence. If the achronal set has no 
edge, this worldline must intersect the set, and it will do it at a point q with 
r{q) — r(g_) = r(g+) > r{q). We have thus proven that for any achronal 
set without edge, the coordinate r cannot be bounded. It follows that any 
achronal set in the manifold cannot be both compact and without edge. 

Finally, let us remark that in Ref.2 there appear two other singularity 
theorems due to Hawking (marked with numbers 3 and 4, pages 271 and 
272). However, the conditions in these theorems are stronger than those of 
the Hawking-Penrose theorem which we have already proven not to hold in 
the solution. There are some other singularity theorems now available, but 
their hypothesis are mere variations of those here studied. In this sense, we 
believe that our study is somehow exhaustive. In the next section we shall 
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give a brief discussion of the meaning of our results. 

4 Discussion 

We have proven that solution (1) is singularity-free and in accordance with 
the main singularity theorems as well. In fact, once the properties of the 
solution are known, it becomes rather obvious that it is free of singulari- 
ties. Because the solution is globally hyperbolic so that there cannot be 
any Cauchy horizon. However, the congruence of the fluid is trivially com- 
plete and through each point in the manifold there is a worldline of the fluid 
congruence. 

It is also clear that, due to the properties of the solution, any possible 
singularity would have to have some extension and therefore it should man- 
ifest itself in the curvature invariants; but they do not. However, to be on 
the safe side, we have preferred to perform a detailed analysis of which par- 
ticular condition in the different singularity theorems failed to be verified 
in the space-time for the following reasons. First, to illustrate how the so- 
lution can avoid the development of singularities. Second, to get a deeper 
insight in the significance and application of the singularity theorems and, in 
this sense, to show explicitly that neither the energy or causality conditions 
are determinant by themselves for the appearance of singularities, not even 
for solutions filled with reasonable matter everywhere and such that all the 
matter is expanding at a given instant of time. 

As is well known, reasonable singularity-free solutions arc a rarity at the 
moment. This is specially true for solutions with cosmological properties. By 
this we do not mean solutions which can describe adequately the real observed 
Universe (which is not the case here) , but rather solutions which theoretically 
must be considered as cosmological since they have matter everywhere and 
cannot be considered as interiors of some vacuum exteriors (for this it would 
be necessary a timelike surface of vanishing pressure) . Therefore, the question 
arises of how many singularity-free solutions there possibly are and which 
particular properties they must obey. 

It can be thought, for example, that the special equation of state (6) 
would have some importance in order to avoid the singularities. This is 
not true, for in Ref.lO it has been shown that the solution here studied 
belongs to a larger family of singularity-free metrics for perfect fiuids with 
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no equation of state whatsoever. However, all the members in the family 
have in common that they are cylindrically symmetric and, in fact, every 
other solution found in [10] without this symmetry contains singularities. 
Thus, cylindrical symmetry could be somehow important for the avoidance 
of singularities. Of course, we must keep this conclusion as a mere hypothesis. 
On the other hand, the converse of this hypothesis is obviously false, and, 
in fact, if we replace cosh{at) for sinh{at) in (1) we obtain another solution 
of Einstein's equations with similar properties and a well-defined cylindrical 
symmetry but having a big-bang singularity in the finite past [1,11]- 

In any case, it is not clear to us why there should not be general solutions 
(with no symmetry) without singularities if they fulfill the main properties 
shown in this paper, that is to say, if they do not have any kind of causal 
trapped set, which reveals itself as an essential assumption in the singularity 
theorems. 
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